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ABSTRACT 

We employ a high-resolution dissipationless N-body simulation of a galaxy cluster 
to investigate the impact of subhalo selection on the resulting velocity distributions. 
Applying a lower limit on the present bound mass of subhalos leads to high subhalo 
velocity dispersions compared to the diffuse dark matter (positive velocity bias) and 
to a considerable deviation from a Gaussian velocity distribution {kurtosis ~ —0.6). 
However, if subhalos are required to exceed a minimal mass before accretion onto 
the host, the velocity bias becomes negligible and the velocity distribution is close to 
Gaussian {kurtosis ~ —0.15). Recently it has been shown that the latter criterion 
results in subhalo samples that agree well with the observed number-density profiles 
of galaxies in clusters. Therefore we argue that the velocity distributions of galaxies in 
clusters are essentially un-biased. The comparison of the galaxy velocity distribution 
and the sound speed, derived from scaling relations of X-ray observations, results in 
an average Mach number of 1.24. Altogether 65% of the galaxies move supersonically 
and 8% have Mach numbers larger than 2 with respect to the intra cluster gas. 
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1 INTRODUCTION 

Velocity distributions in groups and clusters of galaxies can 
be used to determine their dynamical masses. It is commonly 
agreed upon that line of sight distributions similar to Gaus- 
sian reveal relaxed systems (e. g. C hincarini & Rood 1977; 
iHalhdav et alJ l2004l : iLokas et^ahl l2006t) . Non-gaussianity 
is usually associate d with merging or e ven multiple- 
merg i ng events (e.g. ICoU ess fc DunnI 1199^: Cortese et alJ 
l2004l: lAdami et all 120051: iGirardi et al.l 120051) . Obviously. 



only relaxed systems may yield reliable mass estimates. 
However, cold dark matter (CDM) simulations have revealed 
high velocity dispersions of subhalos compared to the dif- 
fuse dark matter, even if relaxed systems are considered 



jGhigna et alJ2000l:IColm et alj2o6alDiemand et al.l2004D . 

In other words, the subhalo populations show a positive 
velocity bias or are hotter compared to the diffuse com- 
pon ent. For dynamica l mass estimates of observed clusters 
(e.g. iLokas et aLlbOOeTl it is important to know whether the 
velocities of galaxies are biased or not. 

In comparing real galaxy clusters with N-body simu- 
lations, subhalos must be associated with galaxies. There 
have been different selection criteria proposed in the lit- 
erature. As discussed below, the strength of the velocity 
bias depends strongly on the subhalo selection. Therefore 
we compare the velocity distributions of two differently se- 
lected subhalo samples derived from the same N-body clus- 
ter. One sample comprises all bound dark matter substruc- 



tures above a certain mass limit at present time (z = 0) . This 
kind of selection, which has been used in earlier investiga- 
tions, leads to a positive velocity bias. However, the spatial 
distribution of these subhalos is less concentrated than the 
underlying dark matter distribution and not in agreemen t 
with observed galaxy distributions. JPiemand et alJl2004ll . 
A second subhalo sample, with distributions similar to the 
observed galaxy distributions, contains only those subhalos, 
which exceeded a certain mass limit before accretion onto 
the host. Using this accretion time subhalo selection crite- 
rion iNagai fc Kravtsovl i2005l) were successful in matching 
the distribution of obse rved cluster ga l axies with the results 
of N-body simulations. IConrov et alJ 12005) use the maxi- 
mal circular velocity at the time of accretion to assign lumi- 
nosities and stellar masses to (sub)halos and achieve excel- 
lent agreement of modelled and observed galaxy clustering 
properties. This agreement suggests that the luminosity of a 
galaxy is related to the depth of halo potential at the epoch 
of high starforming activity, i.e. to its mass or circular ve- 
locity before entering the group or cluster. 

Semi-analytic modeling of galaxy formation combined 
with high resolution dissipationless galaxy cluster simula- 
tions (Soringcl et al. 2001; Gao et al. 2004) also produces 
similar spatial and velocity distributions for the dark matter 
and the galax ies. Evidently th e subhalos are populated with 
galaxies as in lNagai fc Kravt sov (2005) but in a less trans- 
parent, more model dependent way. Similar spatial distribu- 
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Figure 1. Normalised radial number distribution for the total 
amount of dark matter (dashed line) and subhalos belonging to 
the two different samples a and b displayed in blue dotted and red 
solid lines, respectively. The dashed-dotted lines display subsam- 
ples of sample b. Red and blue colour indicates accretion before 
and after z = 0.15, respectively. The region within 0.2r^i^ is ex- 
cluded from accumulation to reduce the influence of numerical 
overmerging. 



tions of galaxies and dark matter were a lso found in hydro- 
dynamic cosmological simulations (see 'Nagai fc Kravtsovl 
[2005; Sonmier-Larscn ct al. 2005; Maccio ct al. 2QQ^. 

The spatial distribution of tracers in cosmological dark 
matter halos is related to their velocity distribution to a 
very good approximation via the sphe rical, stationary Jeans 
equation jPiemand et aljl2004l . l2005h . A spatially extended 
component (like subhalos) is hotter than the dark matter 
whereas more concentrated subsets (like intra cluster hght 
or globular clusters) are colder. Since cluster galaxies tend 
to trace the total mass distribution, one expects little or 
no difference between galaxy and dark matter velocity dis- 
tribution. Clusters of galaxies are formed by gravitational 
collapse, which results in a nearly Gauss ian velocity distri- 
bution for the diff use dark component (iHoeft et al.ll20oj : 
iHansen et aIll2006D . thus the velocity distribution of galax- 
ies should be close to Gaussian as well. In the simula- 
tion and the two subhalo samples are described. In 
discuss the impact of the selection criterion of the subhalo 
samples on the velocity dispersions. Additionally the aver- 
age Mach number of galaxies orbiting in the intra cluster 
medium (ICM) is derived. In i^we summarise our results. 



2 SIMULATION AND SUBHALO SAMPLES 

We analyse a cluster-sized dark matter halo generated 
within a cosmological N-body simulation {^m = 0.268, 
SIa = 0.7, (Tg = 0.7, /iioo = 0.71). The peak mass resolu- 
tion is 2.2 X 10^ M(7i with a softening le ngth of 1.8 kpc. The 
cluster ("D12" in lDiemand et al]l2004l) has a virial mass of 



3.1 X lO^^'M© at z = which corresponds to ~ 14 x 10*^ 
particles within the virial radius of 1.7 Mpc ^ . 

We create two different subhalo samples. On the one 
hand we use SKID^ with a linking length of 5 kpc and iden- 
tify all bound structures comprising at least 10 particles as 
subhalos. This way we find 4239 subhalos within the virial 
radius of the cluster at 2: = 0. Subsequently, this subhalo 
sample is referred to as sample a. On the other hand we 
trace back the most bound particle of the subhalos, which 
were identified by SKID in the same manner as mentioned 
above, and compare their positions with those of field ha- 
los which are found with a friends-of-friends group finder 
(FOF) in earlier outputs using a linking length of 0.2 times 
the mean particle separation. The final sample encompasses 
only those subhalos, which had progenitors (FOF field ha- 
los) containing a minimum of 200 (4.4 x 10^ M©) particles 
at least once during their field-halo phase. In total, 367 sub- 
halos me et this criterion. These sub halo s are assumed t o host 
galaxies. iNaeai fc Kravtsovl ^2005^ and lConrov eEaDlEooi) 
used this approach to assign galaxies to subhalos derived 
from pure dark matter N-body simulations. Subsequently, 
this sample is referred to as the "galaxy sample" or sample 
b. The galaxy sample is subdivided into two sub-samples 
according to accretion times before and after z = 0.15. The 
old and the young sub-samples contain 174 and 193 sub- 
halos, respectively. We do not intend to assign any stellar 
properties to the subhalos, however it is expected that the 
old sample represents on average redder galaxies, since star 
formation in these galaxies may be efficiently suppressed by 
interactions with the dense intra cluster medium. 

Figure displays the cumulative radial number distri- 
bution for the different samples and the diffuse dark mat- 
ter component. The location of the most bound particle is 
chosen as centre, which is assumed to coincide with the cen- 
tral brightest galaxy of observed clusters. We start the cu- 
mulation for all components, subhalos, galaxies and diffuse 
dark matter, at 20% of the virial radius (0.2r„ir)- The in- 
ner 20% is excluded because the high density environment is 
likely to artificially remove substructure by numerical over- 
merging. Moreover the survival of galaxies in this very inner 
region also depends on the mass distribution within their 
inner, baryon dom inated parts llNagai fc Kravtsovl lioO^ : 
iMaccio et al.ll20o3) . The galaxy sample profile (solid line) 
and the diffuse dark matter profile (dashed line) are very 
similar, whereas the halos of sample a (dotted line) show a 
definite deviation. The splitting of the galaxy sample accord- 
ing to accretion times, leads to a strongly concentrated old 
subsample (red dotted-dashed lines) . The young sub-sample 
(blue dotted-dashed lines) is more similar to sample a. 

The selection for sample a based on the current mass 
of the subhalos ignores the history of the individual subha- 
los. Thus a recently accreted low mass halo and a tidally 
striped old (i.e. early accreted) subhalo are treated equiva- 
lently. Due to the steep mass function of field halos, most 
systems ever accreted had masses not much larger than the 
minimal mass. Those halos that who still lie above this min- 
imal mass today are mostly systems which have lost little 



^ According to the definition used here the virial radius encloses 
368 times the mean matter density. 

^ http:/ /www-hpcc. astro. washington.edu/tools/skid. html 
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Figure 2. Normalised 3D velocity dispersion profiles {vc 
958kms~^) for all dark matter particles (dashed line) and halos 
belonging to the subsamples a and b displayed in dotted and solid 
lines, respectively. The error bars display Icr deviations. The re- 
gion within 0.2r„ir is likely to be affected by overmerging. Open 
stars and squares indicate the subdivision of sample a in young 
and old galaxies, respectively. The solid green line shows the ve- 
locity dispersion derived from X-ray temperature profiles for clus- 
ters of comparable size. 



mass, i.e. r ecently accreted halos in the outer part of the 
cluster fsee iKravtsov et aLlliool IZentner et al.ll2nO,'Th . The 

selection criterion for the galaxy sample ensures that only 
subhalos with a substantial initial mass (200 particles) are 
counted as members of the sample. Since these halos must 
have retained at least 10 particles to be found by the SKID 
halo finder at z = 0, they can loose 95% of their initial mass 
on their orbits within the cluster potential well (or even more 
if they were more massive at the moment of accretion). The 
subhalos of the galaxy sample are durable. In that respect 
the galaxy sample is very similar to the diffuse dark matter 
component, which by default is indestructible. 



3 DEPENDENCE OF THE VELOCITY 
DISTRIBUTION ON THE SUBHALO 
SELECTION 

Figure |5| displays radially binned velocity dispersions of the 
samples a and b with open and filled symbols, respectively 
and the diffuse dark matter indicated by the dashed line. All 
dispersions are in units og the circular velocity of the host 
halo Vc,max = 958kms~^. The dotted vertical line marks 
the central region which is prone to numerical overmerging. 
Table^displays the characteristic values of the velocity dis- 
tributions excluding the inner 0.2r„ir. The qualitative pic- 
ture does not change if the central volume is included. All 
velocities are computed relative to the centre of mass ve- 
locity {vcom), the average velocity of all particles within 
rvir- The dispersion of the diffuse dark matter component is 
shown as black dashed line. The solid green line shows the 



(^/^c.max 

kurtosis 

p 

N 



0.665 ±0.007 
-0.61 ± 0.16 
-0.073 ± 0.031 
4115 



0.592 ±0.023 
-0.16 ±0.44 
0.17 ±0.11 

339 



0.5816 ± 0.0001 
-0.241 ±0.003 
0.1103 ± 0.0005 
10.9 X 10^ 



Table 1. Velocity dispersion, kurtosis (v^/cr* — 3) and anisotropy 
parameter (3 = 1 — 0.5(Tjjjjj/(t^^j of the sub-halo sample (a), the 
galaxy sample (b) and the diffuse dark matter (diff). 



velocity dispersion derived from X-ray temperature profiles 
fo r clusters with compar able size using the relati ons given 
in lVikhlinin et al.1 ll2005l) and lEvrard et al.1 lll99d) . 

For all radii the velocity dispersion of sample a devi- 
ates by more than ~ 15% from the diffuse component. The 
difference between these two distributions increases towards 
the centre, approaching values as large as ~ 30% at 0.2rvir- 
There also appears a slight deviation of the velocity disper- 
sion of sample b compared to the diffuse component, how- 
ever, these deviations lie within the la uncertainty range. 
The velocity dispersion profile of the galaxy sample (sample 
b) and the diffuse dark matter component are very similar. 
The old and recently accreted galaxy subsubsamples reveal 
lower and higher velocity dispersions, respectively. The ve- 
locity dispersions of the young subsample agree with the 
total galaxy sample at large radii but show excess towards 
the centre. Dispersions of the old galaxy subsample devi- 
ate below the total galaxy sample at large radii. The mean 
velocity dispersions of the subsamples compared to the ve- 
locity dispersions of all galaxy halos are add = 0.95aaii and 

anew ~ 1.04o-aH. 

As discussed before, the similarity of the galaxy sam- 
ple and the diffuse component in the density profiles can 
be explained by the long lifetime of the sample members. 
Therefore halos of sample b can be considered as a set of 
durable particles similar to the simulation particles, but with 
much larger masses. Due to energy conservation (and with- 
out including dynamical friction or other energy redistribut- 
ing mechanisms) the gravitational collapse of a distribution 
of different mass particles initially leads to equal velocity 
dispersions within different mass bins. In this scenario nei- 
ther spatial nor velocity biases are expected between sam- 
ple b and the diffuse component. On the other hand Fig- 
ure |21 indicates a prominent velocity bias or offset between 
the diffuse component and sample a. The average lifetime 
of sample a is shorter compared to sample b. Sample a is 
weighted toward subhalos that recently entered the host halo 
and consequently move faster. This mechanism shifts the 
average velocity of the remaining subhalos in the sample to- 
wards higher values and causes a positive velocity bias (s ee 
'Ghigna et alfeOOOtlCoh'n et al.l200ol:lDiemand et alJl2004^ . 
Figurc|3|compares the projected velocity distributions of the 
two subhalo samples and the diffuse component. The dot- 
ted, solid and dashed lines are the Gaussian distributions 
derived from the dispersions of the respective components. 
Table displays the actual values of the velocity dispersions 
in units of the maximum circular velocity of the host halo 
(vc,max = 958kms~^). The distributions of the galaxy sam- 
ple and the diffuse dark matter are very similar. The velocity 
dispersions and kurtosis of these two samples agree within 1 
a uncertainty. The velocity dispersion of sample a, however. 
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Figure 3. Velocity distributions of sample a, b are displayed 
with blue open and red filled circles. The diffuse dark matter is 
indicated by black triangles. The lines give the Gaussian distri- 
bution according to the mean velocity dispersion of each sample 
Only objects in the range 0.2 < r/r^ir < 1 are used. 

exceeds the two others by ~ 15%. The corresponding veloc- 
ity distribution (open circles, Figure |^ is flat-topped and 
not well approximated by a Gaussian distribution (dotted 
line). There appears to be a lack of slow moving subhalos 
and a slight excess of high velocity subhalos causing a neg- 
ative kurtosis of ~ —0.6. These features can naturally be 
explained by the loss of earlier accreted, slow moving, sub- 
halos due to tidal truncation. In this context, loss means 
dechne in particle numbers below the resolution limit of 10 
particles. 

The shapes of observed galaxy velocity distributions in 
relaxed clusters can in principle be used to infer if they follow 
a biased, flat-topped or unbiased, Gaussian distribution. In 
practice this are diflicult measurements since large numbers 
of galaxies and careful removal of interlopers are needed to 
achieve a signiflcant result. T here are some first hints for fl at- 
topped velocity distributionsi fvan der Marel et al] JMOO') re- 
port a negative hi = —0.024 ± 0.005 (which is com para- 
ble to the subhalo samples, see lDiemand et 81112004^ after 
stacking 16 CNOCl clusters and excluding the cD galaxies. 
[Lokas ct al. (2006) found negative kurtosis values (around 
—0.4) in 5 out of 6 nearby relaxed Abell clusters. However, 
the deviations from a Gaussian are only about Icr for these 
5 individual systems. Our study suggests that a kurtosis 
which is significantly more negative than the one for the dif- 
fuse component (which has k ~ —0.15) could be an indicator 
for positive velocity bias (and a related spatial anti-bias) . It 
need not necessarily be related t o tangential orbits (nega- 
tive /3) as often ass umed (see e.g. Ivan der Marel et al.ll2000l 
based on models of lGerhardlFlQQ^h . 

The orbital anisotropy of the galaxy sample is not sig- 
nificantly different from the dark matter background, both 
are slightly radial in cluster D12. Sample a on the other 
hand shows marginally tangential orbits. Note that there is 
signiflcant variation from halo to halo in the [3{r) profiles. 



M > 1 A/ > 2 M >3 

64.44 % 8.33 % 0.18 % 

Table 2. The percentage of galaxies which are expected to exceed 
the Mach numbers listed in the upper line. 

The average over six relaxed clusters similar to D12 shows 
that the dark matter /3(r) grows from zero (i.e. isotropic) 
to about 0.35 near the virial radius and total subhalo pop- 
ulations (corresponding to our sample a) show a similar be- 
haviour with a weak tendenc y to be closer to isotropic on 
average dPiemand et al.ll2004h . 

The green line in Figure |5| displays the velocity dis- 
persion (temperature) profile of X-ray gas in clusters 
with masses compara ble to the cluster analyzed here (see 
IVikhlinin et al.l|2005^ . Despite all the complex gas physics 
involved, it is very similar to the diffuse dark matter and 
galaxy sample profiles. This finding can be used to estimate 
the typical Mach numbers of galaxies with respect to the 
ICM. Assuming adiabatic sound speed {vs = \/ ycrl, where 
(Ji is the ID velocity dispersion of the system and 7 = 5/3 is 
the adiabatic constant) and integrating the Gaussian distri- 
bution of the galaxies results in an average Mach number of 
~ 1.24. The distribution of Mach numbers is displayed in Ta- 
ble|5| For supersonic galaxy motions lea ding bow shocks and 
ram pressure stripped tails are present dStevens et al 
which can be detected by X-ray observations. Tails of super- 
sonic galaxies are expected to be more irregular th an these 
of subsonic moving galaxies (iRoedieer et al ] l2006ll . More- 
over, since ram pressure is proport ional to the ICM de nsity 
times the galaxy velocity squared jGunn fc Gottill972h . the 
appearance of leading bow shocks reduces ram pressure and 
decreases the stripping efficiency, which may have impact on 
the abundance profiles in groups and clusters. 



4 SUMMARY AND CONCLUSIONS 

Using a cold dark matter simulation of a cluster sized 
(3.1 X 10^* Mq at 2 = 0) host halo, we find that differ- 
ent subhalo selection criteria change the resulting velocity 
distributions. We analyse the velocity distribution of two 
differently selected subhalo samples. Sample a contains all 
presently found subhalos with masses above 2.2 x 10* Mq 
(10 particles). Similar selection criteria have commonly been 
used for investigations of the velocity bias in N-body sim- 
ulations, however they do very likely not generate subhalo 
samples which are comparable to galaxies in groups and clus- 
ters. Sample b comprises subhalos, which were able to ac- 
cumulate more than 200 particles before entering the host 
halo. This kind of selection results in number density pro- 
files which are simil ar to galaxies observed in groups and 
clusters (see iKravtso v et al. 2004; Nagai & Kravtsov 2005j). 
Our main conclusions are as follows: 

(1) In agreement with other authors, we find an enhance- 
ment of the velocity dispersion in the range from 15% to 30% 
if sample a is compared to the diffuse dark matter compo- 
nent. On average sample a comprises more recently accreted, 
fast moving, sub-halos since the early accreted, somewhat 
more slowly moving, halos are prone to tidal dissolution. 
The positive velocity bias in sample a results from a lack 
of slow moving sub-halos, i.e. a flat topped non Gaussian 
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velocity distribution with negative kurtosis k = —0.6. 

(2) We find no significant velocity bias between sample 
b and the diffuse component. Both have nearly Gaus- 
sian velocity distributions (fc ~ —0.15) and small radial 
anisotropics (/? ~ 0.15). Since sample b resembles the spatial 
distributions of galaxies within clusters, it seems reasonable 
to identify sample b with such galaxies. We conclude, that 
the velocity distribution of cluster galaxies is very similar to 
the underlying dark matter velocity distribution. This find- 
ing supports the assumption of not applying a spatial or 
velocity bias when estimating th e cluster masses via galaxy 
kinematics (see e.g. iLokas et ahl EoOSl. 

(3) The difference between sample a and b lies in the life- 
times of the subhalos. Many subhalos of sample a are low 
mass objects, lying only a little above the mass limit. If 
these subhalos lose a small fraction of their mass due to tidal 
forces, they will no longer be members of the sample. On the 
other hand members of sample b have to lose at least 95% of 
their mass to be removed from the sample. Likewise, massive 
galaxies in clusters are assumed to survive for a long time 
after accretion. This can explain similar properties of the dif- 
fuse dark matter, sample b and luminous galaxies in clusters. 
However, this picture may change in smaller and/or older 
host systems where dynamical friction and tidal forces are 
mor e important. For insta nce, fossil groups are presumably 
old dP'Onghia et alJl2005l) and may have turned a substan- 
tial fraction of their old, sl ow moving, satellite galaxies into 
diffuse intra-group ligh t ijD^Rocha fc Mendes de Qliveiral 
I2OO5I: iFaltenbacher fc Mathews! l2005l) . It is expected, that 
the spatial and velocity distributions of fossil groups show 
similar features (flattened central number-density profile 
and a lack of slow moving galaxies) as found in sample a. 

(4) The mean velocity dispersions of the whole galaxy sam- 
ple compared to the old galaxy subsample differ by 5%, thus 
the resulting mass estimates based on would deviate by 
a factor of 10%. A similar trend is found in observations, 
if more recently accreted galaxy populations are included 
for the computation of the total velocity dispersion (see 
e.g. iMendes de Oliveira et"alll200d) . 

(5) We find an average Mach number of 1.24 for 
galaxies moving within a relaxed cluster (compare to 
IFaltenbacher et alJl2005ll . Altogether 65% of the galaxies 
move supersonically and 8% show Mach numbers larger than 
2. The appearance of shocks affects the interaction between 
galaxy and cluster gas in various ways. In particular shocks 
ahead of supersonic moving galaxies reduces the ram pres- 
sure exerted on the gas in their disks. 
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